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Abstract. 

We investigate the initial- value problem of the non-linear Liouville hierarchy. For the general 
form of the interaction potential we construct an explicit solution in terms of an expansion over 
particle clusters whose evolution is described by the corresponding-order cumulant of evolution 
operators of a system of finitely many particles. For the initial data from the space of integrable 
| functions the existence of a strong solution of the Cauchy problem is proved. 
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1 Introduction 

t-H ! 

The nonlinear Liouville hierarchy that describes the evolution of correlation functions, arises 
in many problems of statistical mechanics concerning many-particle systems [HEJCJ- However, 
today it is still insufficiently studied from the mathematical point of view. 

It is well known [1], that all possible states of a classical system of a finite number of particles 
are described by the functions interpreted as probability density functions. These functions are 
solutions of the initial- value problem of the Liouville hierarchy — the first-order partial differential 
equations, whose characteristic equations are Hamilton equations. If the state of a system is 
presented in terms of a cluster expansion in new (correlation) functions one evidently obtains an 
^ ■ equivalent description of this state. Now evolution of the correlation functions is determined by 

the nonlinear Liouville hierarchy — certain nonlinear first-order partial differential equations. 

In this paper an explicit solution of such nonlinear equations is constructed and presented as an 
expansion in terms of particle clusters whose evolution is described by a cumulant (semi-invariant) 
of the evolution operators. The latter are determined by the solutions of the characteristic 
equations of the linear Liouville equation, i.e., the Hamilton equations. The interaction potential 
of the general form is considered, which makes it possible to describe the general structure of the 
generator of the nonlinear Liouville hierarchy. The existence of a strong solution of the Cauchy 
problem with initial data from the space of integrable functions is proved. A formal treatment 
of the nonlinear Liouville hierarchy for the case of a pairwise interaction potential, was given by 
Bogolyubov and Green [HE]- 

It should be noted that the nonlinear Liouville hierarchy is basic in the substantiation of 
the derivation of the nonlinear Bogolyubov hierarchy [2] whose solutions describe the correlation 
dynamics of infinite systems of particles. Moreover, this concerns the mathematical substantiation 
of the correlation-weakening principle. The correlation functions may be employed to directly 
calculate the specific characteristics of the system, i.e., fluctuations, defined as the average values 
of the square deviations of an observable from its average value, as well as macroscopic values 
which are not averages of observables. The construction of the nonlinear Bogolyubov hierarchy 
and the analysis of the solutions thereof for on the basis of the results obtained here will be given 
in an other paper. 
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2 Initial- Value Problem of Nonlinear Liouville Hierarchy 

Let us consider a system of non-fixed number of identical particles with phase coordinates X{ = 
iliiPi) £ xR", > 1. The relevant Hamiltonian is given by the formula 

n 2 n n 

#n=E^"+E E *fc(9ii,9i 2 ,---,9i fc ), 

i=l fc=l ii<J2<...<ifc=l 

where $fc is the fc-th-order interaction potential. In what follows we assume that the interaction 
potential k > 1, satisfies the necessary conditions which provide the existence of global in 
time solutions of the Hamilton equations. 

The state of the system can be described by a sequence g(t) = (0, gi(t,x±), . . . 
. . . , g n (t, xi, . . . , x n ), . . .), of correlation functions g n (t,xi, . . . ,x n ) defined on the phase space 
W n x R^ n , v > 1 and symmetric with respect to the permutations of arguments x±, . . . ,x n . 

The evolution of the states of the above system is described by the initial-value problem of 
the nonlinear Liouville hierarchy 



j t g n (t,Y) = (-C lYl (Y))g n (t,Y) + 



|p| 

+ E E (-^ (U z 0) n (!) 

p ; y=MXiZ l CX i I U z i\ i=l XiCP 

IPI>1 

5n(t,n| t=0 = 5n(0,y), n>l, (2) 

where the following notation is used Y = (x\, . . . ,x n ), \Y\ = n denotes the number of elements 
of the set Y, Y 1S the sum over all possible (in this case) decompositions P of the set Y into 
P:V U v. 

i 

|P| nonempty mutually disjoint subsets, Y ' ^ s the sum over all nonempty subsets Z; L C X; L . 
The Liouville operator C n for the Hamiltonian H n is described by the formulas 

n q n n 

C"n = £n{ x li • • • > %n) = ^ ] {Pi> ~K ) ^ ] ^ ] ^fe (^ii j ^2 J ■ ■ ■ > x ik)i (^) 

i=l yj fc=2 ii<i 2 <...<ifc=l 

C k (Xjj , Xj 2 , . . . , Xj fc ) = — y~^(g^~^fc(gii) gi 2 > • • • 5 9ifc)> ^7^' 

where (■, •) - is the scalar product. 

The simplest examples of the nonlinear Liouville hierarchy are given by: 

^5i(Mi) = -£i(zi)si(*,£i), 

^92(t,xi,x 2 ) = -C 2 (xi, x 2 )g 2 (t, xi,x 2 ) - C l 2 nt (x 1 ,x 2 )g 1 (t,x 1 )g 1 (t,x 2 ), 

—g 3 (t,x 1 ,x 2 ,x 3 ) = -C 3 (x 1 ,x 2 ,x 3 )g 3 (t,x 1 ,x 2 ,x 3 ) + 

+ ( - 4"*(xi, x 2 ) - 4"*(xi, x 3 ) - C 3 nt ( Xl , x 2 , x 3 )) gi (t, Xl )g 2 (t, x 2 ,x 3 ) + 
+ ( - 4">i,x 2 ) - 4">2, x 3 ) - £ 3 n \x u x 2 , x 3 )) gi (t, x 2 )g 2 (t, xi,x 3 ) + 
+ ( - 4"*(xi,x 3 ) - C 2 n \x 2l x 3 ) - 4"*(xi,x 2 , z 3 ))si(t, z 3 )02(*, zi,x 2 ) + 
-£ 3 n '(xi, x 2 , x 3 )gi (t, xi)g!(t, x 2 )gi(t, x 3 ). 
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We note that in the case of a pairwise interaction potential, (k = 2), the nonlinear Liouville 
hierarchy ([TJ) is simpler. For example, expression for g 3 {t) does not contain members with C 3 nt . 
This case was considered by Green [2J. 

In this work the solution (p])-© is shown to be given by the formula 

g n (t,Y)= ^1^^)11^1(0,^), (4) 

i 

where Y = (xi, . . . ,x n ), Yx t = (X\, . . . , X\ P \), Yl ~ the sum of all possible (in this case) 

P:Y \JXi 

i 

decomposition P of the set Y into |P| nonempty mutually disjoint subsets Xi. The evolution 
operator 2tipi(i) i.e., the cumulant (semi-invariant) of the order |P| is given by the expression 

*\*\(t,Y Xi ) = E (-l) |P ' hl (|P'|-l)! II S\z H \(-t,Z k ). (5) 

P'-Y x .=\JZ k Z k c?' 

k 

The evolution operators S n (—t), re > 1, are given by 

(S n (-t)g n {0))(x 1 , ...,x n ) = g n (0,X 1 (-t,xi, . . . ,x„), . . . ,X n (-t,x 1 , . . .,x n )), (6) 

where {Xi(—t, x%, . . . , x n )}f =1 is the solution of the relevant initial- value problem for the Hamilton 
equations. The properties of a group of evolution operator ([6]) are described in |6j. 

Let us consider the simplest examples of expansions (Jl|) with the following notation: the 
argument x\ U Xj implies, that two particles i-th and j-th evolve as a cluster. Thus, if the 
arguments of the operator are clusters, they enter on equal terms the expansions in series of the 
evolution operators S n (—t), the order of the cumulant being equal to the number of its cluster 
arguments . 

gi(t,xt) = %i{t,xi)gi(Q,xi), 

g 2 (t,xi,x 2 ) = 2ti(t,xi U x 2 )g 2 (0,x 1 ,x 2 ) + 2l 2 (i, %i, x 2 )gi(0, xi)5i(0, x 2 ), 

g 3 (t,x 1 ,x 2 ,x 3 ) = QLi(t,xi U x 2 U x 3 )g 3 (0,xt,x 2 ,x 3 ) + 2t 2 (£, sci, x 2 U £3)51(0, £i)5 2 (0, x 2 , x 3 ) + 
+% 2 {t,xi U x 3 ,x 2 )gi(0,X2)g 2 (p,x 1 ,x 3 ) + % 2 (t,x 1 U x 2 , x 3 )gi (Q,x 3 )g 2 (0, x 1: x 2 ) + 
+2l 3 (t, x 1 ,x 2 ,x 3 )g 1 (0, 2:1)51(0, ac 2 )5i(0, x 3 ), 

where, for example, the cumulants ((5|) are given by 

2l 2 (i, xi U x 2 , x 3 ) = S 3 (-t, xx, x 2 , x 3 ) - Si(-t, x 3 )S 2 (-t, xi,x 2 ) 
% 3 (t,xi,x 2 ,x 3 ) = S 3 (-t,xi,x 2 ,x 3 ) - 

-Si(-t,xi)S 2 (-t,x 2 ,x 3 ) - Si(-t, x 2 )S 2 (-t, x±, x 3 ) - Si(-t, x 3 )S 2 (-t, xi, x 2 ) + 

+2!5i(-t, xi)Si(-i, x 2 )5i(-t, x 3 ). 

Formally, the nonlinear Liouville hierarchy ([1]) can be derived from the sequence of (linear) 
Liouville equations which describe the evolution of all possible states of the system of non-fixed 
number of particles (the sequence D(t) = (l,Di(xi), . . . ,D n (xi, . . . ,x n ), . . .), with the function 
D n (t) being regarded as the density of probability distribution of the re-particles system) [I] 

^1 --CD (t) 
D n (t)\ t=Q =D n (0), n>l, 
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provided the state of the system is described in terms of correlation functions, i.e., 

9n (t,Y)= Y, (-l)' Phl (|P|-l)! II D \Xi\(fM (7) 

?:Y=\JXi X,CP 

i 

or at the initial time instant 

g n (0,Y) = Y (-l) |Phl (|P|-l)! II DmfrXi), n>l, (8) 

?:Y=\JXi XiCP 

i 

where Y = (xi, . . . ,x n ), Yl i s the sum over all decompositions P of the set Y into |P| nonempty 

p 

mutually disjoint subsets Xi. 
For example, 

gi{t,xi) = D 1 (t,xi), 

92(t,xi,x 2 ) = D 2 (t,x 1 ,x 2 ) - D 1 (t,xi)D 1 (t,x 2 ). 

The solution ^ of the initial-value problem (JXJ)- (EJ) can be formally derived from© and ([ED 
provided one takes into account that, within the context of (j5|), we have 

D lXil (0,X i )= Y II 9\z k \(0,Z k ). 

P2:Xi=\JZ k Z k C?2 

k 

Then, inasmuch as the solution of the Liouville equation is given by 

D n (t) = S n (-t)D n (0), 

where S n (—t) is determined by the solutions of Hamiltonian equations according to (|6|), we have 
9n (t,Y)= Y, (-l)' Plhl (|Pl|-l)! II X 

i 

xS\ Xi \(-t,Xi) Y II 9\z k \(P,Z k ). (9) 

k 

Having collected in j9|) the terms with similar product of functions g\z k \(0, Zk), one obtains (jl]). 

3 The existence theorem for the initial- value problem of the nonlinear 
Liouville hierarchy. 

Suppose L\ is the Banach space of integrable functions g n (xi, ■ ■ ■ , x n ), defined on the phase 
space x W 11 , v > 1 of n-particle system, symmetric under the perturbations of arguments. 
The norm of an element g n of L\ is denoted by 

\\g n \\ = / dxi . . . dx n \g n (xi, . . . ,x n )\, 

L\ C L\ is a subspace of continuously differentiable functions with compact supports. 
The following theorem is true. 
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Theorem 1. // g n (0) € L\q C L\, n > 1, i/ien /or iGl 1 i/iere earisis a unique strong solution 
to the initial-value problem (OP-0 of the nonlinear Liouville hierarchy J7P given by 

g n (t,Y) = %\P\& Y Xi) II WQM* 

P:Y=\JX t X t CP 

i 

where 2tipi(i) is a cumulant (semi-invariant) of the order \P\ of the evolution operators 

flonasamejibcmeo. Let us show that the expansion ([4]) is defined in L\. 

Indeed, within the context of the corollary of the Liouville theorem [BJ (isometric property of 
operators S\z k \(—t, Z k )) the following estimate is valid 



\9n 



f dY\ Y, II 

P:Y=\JX, X,CP 



= l dY \ e e (-i) ipVi (ip'i-i)! n %z k \{-t,z k ) n w°>*oi < 

P:Y=[jX lP '..Y Xt =[JZ k Z k CP' ^CP 

1 fc 

< E E (i p 'i - !) ! II n»(°)ii ^ n!e " +1 E II n»(°)ii < °°. 

p : y=U^p' : y x .=uz fc ^cp p : y=u^^cp 

i.e., #n(£) € for any i £ M 1 . 

Let us prove that the expansion (|1J is a strong solution of the Cauchy problem of the nonlinear 
Liouville hierarchy (EIJ) - d2]) - 

To do this we first differentiate the functions g n if) with respect to time with regard for the 
point-by-point convergence. Let g n (0) € L^ , then inasmuch 

jj.S\Xi\(—ti Xi) = —£\Xi\(Xi)S\Xi\(—t,Xi), 

where the Liouvilian C\xa is defined by (EI), and according to J9j), for each fixed point Y on any 
compact from W n x R un , we have 

| 5n (t,Y) = | a |P| II 9\x 1 \(0,X l ) = 

P:Y=\JXi XidPi 

i 

=i e (-i) |piM (iPii-i)! n e n 9\ Zk \(o,z k )= 

PV.Y=\JX, X.CPi P 2 :X 1 ={JZ k Z k CP 2 

i k 

= y (-i) |Pihl (iPii-i)! E (-*w*i)) n w-«)x 

P v .Y=\JXi X,CPi XiCPi 

i 

X E II 9\z k \(0,Z k ). 

P 2 :X 1 =\JZ k Z k CP 2 

k 

In view of the equality 

S\ x .\(-t,Xi) Y U 9\z k \(0,Z k )= Y II 9\z k \&Z k ), 

P 2 :Xi=\JZ k Z k CP2 P 2 :Xi=\JZ k Z k CP 2 

k i 



5 



which follows from (El), we have 



d 



j t 9n(t,Y)= e (-i) |pihi (iPii-i)! e n e II W*>^) 



= (-£| Y |(y)) 5n (y) + E E (-i) |P2hl (|p 2 | - 1)! x 

fvY={jXiP 2 :Y x .=\JZ k 

i k 

Pi>i 

x E ("% fc |(^)) II aix^Xi), 

Z k CP 2 XiCP! 

the notation being similar to that in and Yx t = • • • ,^|p|) is the set whose elements are 
|P| subsets Xi CY. 

Since the following identity holds 

E (-i) |phl (|p| - 1)! E = E ( - ^ ( U 



*CP 



U ^1 i=i 



where is a sum over all nonempty subsets Zi C X,- L and operator is denned by 



Zi<ZX{ 



'•int 

I U Zi\ 



formulas we come to equations ([T]). Thus, formula (Jl]) determines the solution of the initial- 
value problem of the nonlinear Liouville hierarchy JIJ) from the viewpoint of the point-by-point 
convergence. 

Let us show that the strong derivative of the solution fll}, reproduces the generator of the 
nonlinear Liouville hierarchy ([T]) in the subspace C L\ . 

For <7n(0) E L\ , in the sense of convergence norm of the space , we have 

E E (-i/'-Wi-i)^ 1 

P:Y=\JX iP '. Yx .=VZ k 



lim 

At^O 



At 



x [ n s lZkl (-(t+At),z k ) n six^o,^)- n n ^ko,^) 



Z k CP' 



X,CP 



Z k CP' 



X,CP 



e {-crnfr)) n n =°- 



ZiCP 



XiCP 



Indeed, within the context of the corollary of the Liouville theorem and the group property of 
operators S n (— i), n > 1, ^ we obtain 



lim 

At^O 



E E (-i) |phl (|p'i-i)! 

P:Y=\JX iP > :Yx .=UZ k 



At 



lim 

At^o 



X,CP 



[ n s\z k \(-(t+At),z k ) n 5 | X .|(o,Xi)- n s lZkl (-t,z k ) n 

Z k CP' XiCP Zk CP' 

e i-c\ Zi \{Zi)) n s\ Z .\(-t,Zi) n 51^1(0,^)] i = 

ZiCP 1 ZiCP' X »CP 

E E (-1) |P ' | - 1 (|P'|-D![^[ II S lZkl (-At,Z, 

P:Y=\JX iP '., Yx .=]JZ k Z k CP' 



k) X 



x n ^,(0,^0- n ^i(o,^)j - e (-4^1(3)) n w^,* 

XiCP XiCP ZiCP' 
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H E E (-d 

p : y=U^p':Y x,=UZ fc 



IP -1 



(IP 



1)! lim 



[ill n 



Z k CP' 



x n 5 | Xi |(o,x 4 )- n 9^0, Xi)\ - y i-twfr)) n ^(o,^ 

^iCP XiCP ZiCP' x ^ p 

The last equality is valid becouse the integrand in this expression tends to zero as At — > 
uniformly in Y on any compact set. Therefore, we can pass to the limit as At — > in the 
integral. Hence, equality (|IJ is differentiable in the norm of the space L\. Thus, for the initial 
data from C L\ the corresponding Cauchy problem has a unique strong solution that given 
by expansions CD)-©- □ 



4 Properties of the solution of the nonlinear Liouville hierarchy. 

We introduce the evolution operator of the solution Q 

E a ip|(*> y *i) II wQM = (a*0/(o))) m (y). (io) 

P:Y=\JXi XiCP 

i 

For the evolution operator (fTOj) . the group property is valid, i.e., 

(a tl (2i t2 (<?(()))) ) |y| (Y) = (2t t2 (2t tl ( ff (o)))) |y| (y) = (2i 4l+42 ( 5 (o))) |y| (y). 

Indeed, for g n (0) € , n > 1 and for any ii, *2 G K 1 , according to the notation Q and JS]), we 
have 

(2t tl (2t t2 ( 5 (0)))) |y| (y) = 

= y »iPi(*i. y jfi) n e %i(*2,(^h) n ^i(°.^)= 



« i 

= e e (-i) |piM (iPii-i)! n *w-*i.<w n >< 

PiVHJ^Pi^jr^UQfc Q fc CP' ^CP 

x e e (-i) |p2M (ip 2 i-i)! n s lRkl (- t2 ,R k ) n mi {o,z{). 

p':X l =\JZ l ^2.{X i ) z =\JR J R k CP 2 Z,CP' 

! 3 



Having collected the items at identical products of the initial data <? n (0), n > 1, and taking into 
account the group property of the evolution operators S n (—t), n > 1 ([6]), we obtain 

(2l il (2l t2 ( 5 (0)))) |y| (y) = 

= E E (-1) |P ' M (|P'|-1)! II %|(-(ti+* 2 ),^) n 9\x t \(0,X l ) = 

P.Y=\JX lP ':Y x .=[JZ l Z t CP' ^(ZP 

i 1 I 

= Y %\*\<ti + t2,Y Xi ) n g lXi \(0,Xi) = (a tl+t2 ( 5 (o))) n (y). 

P:Y={JX t XiCP 

i 

Similarly, 

(2i t2 (2i tl ( 5 (o)))) n (y) = (2i tl+t2 ( 5 (o))) n (y). 
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Let us consider the property of the solution ^ for one physically motivated example of 
the initial data that is to say if the initial data for Cauchy problem ([XJ)- (EJ) , satisfy the 'chaos' 
condition [3], in other words, the sequences of correlation functions have the form 

s(O) = (l,0i(O,xi),O,O,...), (11) 

Indeed, in terms of the sequences D (0), this condition implies [I] that 

D(0) = (l,D 1 (Q,x 1 ),D 1 (Q,x 1 )D 1 (0,x 2 ),...), 

which implies that particle distributions are statistically independent at the initial time instant. 
Making use of the relation (JHJ), we obtain the initial condition (fTT|) for the correlation functions. 

For the initial data (fTTj) . the formula for the solution (|IJ of the initial- value problem (Op-© 
is simplified and reduces to 

n 

9n(t, X\ , . . • , x n , X\ , • • • , x n 

)l[gi(a,xi)- (12) 

i=l 

In this case the following estimate holds 

\\gn(t)\\<n\e n+1 \\ gi (0)p. 

In the case of the initial data (fTTj) . the solution (fT2]) of the Cauchy problem for the nonlinear 
Liouville hierarchy may be rewritten in a different form. If n = 1, we have 

gi(t,xi) = 2li(t,xi)5i(0,xi) = g\(0,qi -pit,pi) 

i.e., an explicit expression that describes the evolution of one particle. Then, within the context 
of the definition of the 1-st order cumulant, 2li(t), and inverse to it evolution operator 2ti(— t), 
we can express the correlation functions g n (t), n > 2, in terms of the one-particle correlation 
function g\ (t) making use of formula (fT2l . 
Finally, formula (fT2l) for n > 2 is given by 

n 

g n (t,xi, ...,x n ) = %i{t,xi, . . . ,x n ) \\gi(t,Xi), 

i=i 

where %n(t, x\, . . . ,x n ) is the n-th order cumulant of the scattering operators St, i.e., 

n 

St (%1 , ■ ■ ■ j %n) = S n (—t, X\ , . . . , Xn) J^J S\(t,Xi), 

i=l 

whose generators are determined by the operator £ m *([3|) in terms of the interaction potential of 
the system. 
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